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$H ' Abstract 

Corollary 2, Entry 9, Chapter 4 of Ramanujan's first notebook claims that 
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O 



^ • nk \k\ ' 



as a; — )■ oo. This is known to be correct for the case n = 1, but incorrect for n > 3. We show 
that the result is correct for n = 2. We also consider the order of the error term, and discuss 
. a different, correct generalisation of the case n — \. 
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> ' 1 Introduction 

ly^ i Much of Ramanujan's work was not pubhshed during his lifetime, but was summarised in his 

' Notebooks. These were printed in facsimile in 1957 [14], and edited editions have been published 

■ by Berndt [2, 3, 4]. 

, Many of Ramanujan's resuhs were obtained in a formal manner, and he did not state suffi- 

cient conditions for their validity. For example, Example 3, Entry 9, Chapter 4 (page 97 of [2]) 
is 



but clearly some conditions on the function (p{z) are required. Berndt has given sufficient 
conditions, but they do not always hold for Ramanujan's applications of (1). In fact, Ramanujan 
does not claim exact equality in (1), but writes that the left side is "nearly" equal to the right 
side. Thus, some sort of approximation or asymptotic equality is intended. 

To illustrate the use of (1), take (p{z) = x^/r(z+l), where a; > is real. Proceeding formally, 
we obtain 

where 7 = — r'(l) is Euler's constant. 

If equality in (2) is interpreted as asymptotic equality (usually denoted by " ~ ") as x ^ 00, 
then the result is correct. In fact, a classical result (also given on page 167 of [3]) is 

°° (_-\\k—l rfk poo p—t 

^~kl\ " ^""^"^^i ^—dt = 0{e-yx) . (3) 
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Ramanujan's Corollary 2, Entry 9, Chapter 4 (page 98 of [2]) is that, for positive integer n, 
In view of (3), this is equivalent to 




~>-- + 7 (4) 

It is plausible that Ramanujan derived this result from (1) in the same formal manner that we 
derived (2) above, but taking ip{z) = {x^/T{z + 1))". 

Berndt [2] shows that (4) is false for n > 3; in fact, the function defined by the left side of (4) 
changes sign infinitely often, and grows exponentially large as x —> oo. However, Berndt leaves 
the case n = 2 open. 

The aim of this paper is to show that (4) is true in the case n = 2. Theorem 1 in Section 3 
gives an exact expression for the error in (4) as an integral involving the Bessel function Jo{x), 
and Corollary 1 deduces an asymptotic expansion. The most significant term is 0(x~^/^) as 
X — ^ oo. 

Theorem 1 is a special case^ of a formula given on page 48 of Luke [9]. However, the 
connection with Ramanujan docs not seem to have been noticed before. 

In Corollary 2, Entry 2, Chapter 3 of his first Notebook, Ramanujan shows that the function 
on the left side of (2) can be written as 

k 

A;=0 

where = Z^j=i 1/i is a harmonic number {Hq = 0). Thus 

°o ,,,fc / oo r^k 



k=Q ' k=0 



In Section 4 we indicate how Ramanujan might have generalised (5) in much the same way 
that he attempted to generalise (2). 



2 Notation and Preliminary Results 

In this section we give some preliminary results on integrals involving Jq. These results may be 
found in the literature, but for completeness we sketch their proofs. 
Recall that 

k=0 

Jo{z) is an entire function, but we are only concerned with its behaviour on the positive real 
axis. For small positive x, Jo{x) = 1 + O(x^). For large positive x, Hankel's asymptotic 
expansion [11, 15] gives Jo{x) = 0(x^^/^). These observations are sufficient to show that the 
integral occurring in Lemma 2 below is absolutely convergent. 

In the proof of the following Lemma, (a)fc = a{a + 1) ■ ■ ■ {a + k — 1) = T{a + k)/r{a), and 

F(a 6.c-z) = V^^M^- 
na,b,c,z) ^^^^ ^, 

is a hypergeometric function. 

^Also given in formula 11.1.20, Chapter 11 of Abramowitz and Stegun [1] (the chapter was written by Luke). 
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Lemma 1 For < < |, 



io ' - r(i - M/2) 



Proof. The integral^ is known as Weber's infinite integral [16]. We sketch a proof in the case 
/i < ^, which is all that is needed below. 

Let a be a complex variable. We first evaluate 

1(a) = / e-'^H''-^Joit)dt 
Jo 

for R{a) > 1. Integrating term by term, the power series (6) gives 

7(a) = a-'^r(M)F(|,^;l;-a-2) . (7) 

By analytic continuation, the result (7) holds^ in the right half-plane R{a) > 0. Using a well- 
known hypergeometric function identity^, we deduce that 

m ^ r(l)F(f,f;i;-a^) , a^(-i)F(^,/^;i;-«^)^^ 



r(i-f) V r(^) r(f) J ' 

Since Jo{t) = 0(t~^/^) as t — > oo, our assumption that A* < ^ makes it is easy to justify changing 
the order of integration and taking the limit as a — >■ 0-|-. Thus 

rt^''^Joit)dt= Ihn /(a) = ^^^^Hl) . 

io ^ ^ "^0+ ^ ^ r(i - ^)r{t^) 

The Lemma now follows from the duplication formula for the Gamma function. □ 
Lemma 2 



/ 

Jo 



^"''-^°<^l)* = 0. (8) 



V t 



Proof. A slightly more general result is given in equation 6.622.1 of Gradshteyn and Ryzhik [7], 
and attributed to Nielsen [10]. We show that (8) follows easily from Lemma 1. 
Let n he a small positive parameter. Prom Lemma 1, 



(9) 



Since r(l) = 1, r'(l) = -7, and /xr(/x) = r(l + ji), the right side of (9) is 

The result follows on letting ^ — ?• 0-|-. □ 



^More precisely, the generalisation with Jo(t) replaced by Jv{i): see Sec. 13.24 of Watson [15]. 

generalisation of (7) is given in Sec. 13.2 of Watson [15], and is attributed to Hankel [8] and Gegenbauer [6]. 



See equation (10.16), Chapter 5 of Olver [11]. 
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3 Ramanujan's Corollary for n = 2 

Our main result is the following Theorem, which proves that (4) is valid for n = 2. 
Theorem 1 Let 

OO / 

2k \ k\ 



k=l 

Then, for real positive x, 



e(x) = r '-^dt . 

J2x t 



Proof. Proceeding as on page 99 of [2] , using the fact^ that 
we have 

e{x) = / ^-^^dt / dt + / dt 

^ ' Jo t Ji t Jo t Ji t 

^ Ut+ — dt . 



10 t Jx t 

Now, from Lemma 2 with a change of variable. 



Jo t Jx 



SO 

, , /-"^ Jo(2t) - e"* , /"^ e"* , /■°° Jo(2t) , 
e{x) = / -^^^^^^ dt + / — = / -^^^dt , 

^a; t J X t Jx t 

and Theorem 1 follows by a change of variable. □ 

Corollary 1 Let e{x) be as in Theorem 1. Then, for large positive x, e{x) has an asymptotic 
expansion whose leading terms are given by 

7t\ 13sin(2x + f) ^ ^^^_2^^ 



e(a;) = — ^ ,„ „ ,„ cos [ 2x + — \ -\ — h 0(x 



Proof. Using integration by parts and the fact that xJq{x) + Jq{x) + xJo{x) = (a special case 
of Bessel's differential equation) , it is easy to deduce from Theorem 1 that 

Continuing in the same way, we obtain an asymptotic expansion 

4{2x)^ ,k\k\ Jo{2x)^ ,,fe fc!(fc + l)! 

k=0 k=0 
The Corollary follows from Hankel's asymptotic expansions for Jo{z) and Jq{z) = —Ji{z). □ 



^See page 103 of [2]. 
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4 A Different Generalisation 



An obvious generalisation of (5) is 

as — 7- oo. 

It is easy to show that (10) is vahd for all positive integer n. An essential difference be- 
tween (4) and (10) is that there is a large amount of cancellation between terms on the left side 
of (4), but there is no cancellation in the numerator and denominator on the left side of (10). 
The function {x'^/kl)'^ acts as a smoothing kernel with a peak at ~ x — 1/2. Since 

Hk = lnk + j + Oil/k) , 

the result (10) is not surprising. What may be surprising is the speed of convergence. Brent 
and McMillan [5] show that^ 

oo / oo /r'^N" 



as X ^ oo, where 



\2n 



if n = 1; 

^ " 2nsin2(7r/n), if n > 2. 



In the case n = 2, (11) has error 0(e~^^). Brent and McMillan used this case with x ~ 17,400 
to compute 7 to more than 30,000 decimal places. Prom Corollary 1, the same value of x in (4) 
would give less than 8-dccimal place accuracy''. 

The case n = 3 of (11) is interesting because maxc„ = C3 = 4.5. However, no one seems to 
have used n > 2 in a serious computation of 7. 

Although Ramanujan [3, 13] gave many rapidly-convergent series and other good approxima- 
tions for TT, he does not seem to have given series which are particularly useful for approximat- 
ing 7, except for (3) and (5) above. In his paper [12] on series for 7, he gives several interesting 
series, of which the simplest^ is 

^=1-2 ^(2fe+i) 



- (A; + l)(2fc + l) ' 

but these scries all involve the Riemann zeta function or related functions, so they are not very 
convenient for computational purposes. 

Our analysis has assumed that n in (4) and (11) is a positive integer. It would be interesting 
to consider the behaviour of the functions occurring in these equations for positive but non- 
integral values of n, especially in the range 1 < n < 2. 



References 

[1] M. Abramowitz and I. A. Stegun (editors), Handbook of Mathematical Functions with For- 
mulas, Graphs, and Mathematical Tables, National Bureau of Standards, Washington, 1964 
(reprinted by Dover, 1965). 

We note an error on page 310 of [5]: in the definition of Vp{z), "z/k\" should be "«*/A;!". 
'^More than 15,000 decimal places would have to be used in the computation to compensate for cancellation 
of terms n±{e'^'' /x^) in (4) ! 
*Due to Glaisher: see [12]. 



5 



[2] B. C. Berndt, Ramanujan's Notebooks, Part I, Springer- Verlag, New York, 1985. 
[3] B. C. Berndt, Ramanujan's Notebooks, Part II, Springer- Verlag, New York, 1989. 

[4] B. C. Berndt, Ramanujan's Notebooks, Part III, Springer- Verlag, New York, 1991. 

[5] R. P. Brent and E. M. McMillan, "Some new algorithms for high-precision computation of 
Euler's constant", Mathematics of Computation 34 (1980), 305-312. 

[6] L. Gegenbauer, "Uber einige bestimmte Integrale (June 17, 1875)" , Wiener Sitzungsberichte 
72(2)(1876), 343-354. 

[7] I. S. Gradshteyn and I. M. Ryzhik, Tables of Integrals, Series, and Products, fourth edition 
(trans. Alan Jeffrey), Academic Press, New York, 1965. 

[8] H. Hankel, "Bestimmte Integrale mit Cylinderfunctionen" , Math. Ann. 8 (1875), 453-470. 

[9] Y. L. Luke, Integrals of Bessel Functions, McGraw-Hill, New York, 1962. 

[10] N. Nielsen, Theorie des Integrallogarithmus und verwandter Transzendenten, Teubner, 
Leipzig, 1906 (reprinted by Chelsea, New York, 1965). 

[11] F. W. J. Olver, Asymptotics and Special Functions, Academic Press, New York, 1974. 

[12] S. Ramanujan, "A series for Euler's constant 7", Messenger of Mathematics 46 (1917), 
73-80 (reprinted in [13]). 

[13] S. Ramanujan, Collected Papers of Srinivasa Ramanujan (edited by G. H. Hardy, P. V. 
Seshu Aiyar and B. M. Wilson), Cambridge University Press, Cambridge, 1927. 

[14] S. Ramanujan, Notebooks, two volumes, Tata Institute of Fundamental Research, Bombay, 
1957. 

[15] G. N. Watson, A Treatise on the Theory of Bessel Functions, second edition, Cambridge 
University Press, Cambridge, 1944 (reprinted 1966). 

[16] H. Weber, "Uber einige bestimmte Integrale (Jan. 1868)", Journal fiir Math. 69 (1868), 
222-237. 



6 



